Orbital mechanisms of electron spin manipulation by an electric 

field 

E. I. Rashba 1 !*] and Al. L. Efros 2 
1 Department of Physics, SUNY at Buffalo, Buffalo, New York 14260, USA 
2 Naval Research Laboratory, Washington, DC 20375, USA 
(Dated: February 2, 2008. SpinManipRevSubm.tex) 

Abstract 

A theory of spin manipulation of quasi-two-dimensional (2D) electrons by a time-dependent gate 
voltage applied to a quantum well is developed. The Dresselhaus and Rashba spin-orbit coupling 
mechanisms are shown to be rather efficient for this purpose. The spin response to a perpendicular- 
to-plane electric field is due to a deviation from the strict 2D limit and is controlled by the ratios of 
the spin, cyclotron and confinement frequencies. The dependence of this response on the magnetic 
field direction is indicative of the strenghts of the competing spin-orbit coupling mechanisms. 
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Manipulating electron spins is one of the central problems of the growing field of semicon- 
ductor spintronics ]J that is of critical importance for quantum computing and information 
processing P]. Most of the schemes proposed for computing with electron spins in quantum 
dots (QDs) and quantum wells (QWs) are based on using time- dependent magnetic fields. 
However for applications using time-dependent electric fields instead of magnetic ones would 
be highly preferable, and various mechanisms of spin-orbit interaction l3[ open attractive 
possibilities for electrical control of electron spins. Recently Kato et al. |4( successfully ma- 
nipulated 2D electron spins by a gigahertz electric field. They used a parabolic Al x Ga!_ x As 
QW formed by varying Al-content x = x(z) gradually across the well. The structure was 
specially engineered to achieve gate- voltage control of the g-factor through its dependence on 
x P). Similar data on the electrical control of the g-factor were reported for GaAs/AlGaAs 
(f| and Si/SiGe 7| heterostructures. These achievements pave the way for manipulating 
electron spins in QDs individually. 

Experimental success in achieving dynamical electric manipulation of electron spins raises 
the question about the dominant physical mechanisms controlling the coupling of spins to 
the electric field. The g-tensor modulation resonance technique |4| is based on the different 
dependence of the various g-tensor components on the gate voltage and works when the 
external magnetic field B is tilted to the QW plane. Indeed, under these conditions the 
operators (ergB) and (erg'B), where cr is the Pauli matrix vector and g' = dg/dV is the 
derivative of g with respect to the gate voltage, do not commute. As a result, a time- 
dependent gate voltage V(t) = Vq + V sm(ut + 0) leads to spin flip transitions at the spin 
resonance frequency uj = uj s = gfisB/h, fiB being the Bohr magneton. The concept of 
the g-tensor mechanism of spin-flip transitions suggested using the region of a very small 
g-tensor, \g\ < 0.1, where g is strongly anisotropic and gate- voltage dependent j^f. 

In this letter we develop the theory for a different mechanism of gate- voltage induced spin 
resonance based on the electron's orbital motion and the standard mechanisms of spin-orbit 
coupling. This theory also requires a tilted magnetic field but does not require the g-tensor 
to be small. Just the opposite, the large g factors typical of narrow-gap A 3 B 5 compounds 
are advantageous. Electric dipole spin resonance (EDSR) [3] is especially strong when it 
is excited by an electric field lying in the QW plane. However, we show that it is also 
strong enough in the geometry when the time-dependent potential is applied to the gate, 
i.e., the time-dependent electric field is perpendicular to the well. This geometry is the most 
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suitable for practical devices. Our results demonstrate convincingly that efficient electrical 
spin manipulation can be achieved through the orbital mechanisms of spin- coupling to the 
electric field. 

Two basic mechanisms of the spin-orbit coupling of 2D electrons are directly related to 
the symmetry properties of QWs. They stem from the structure inversion asymmetry (SIA) 



mechanism described by the Rashba termj8 
anism described by the Dresselhaus term 
same order of magnitude 



and the bulk inversion asymmetry (BIA) mech- 
lCj . In GaAs QWs both terms are usually of the 



~y\ while in narrow-gap compounds like InAs the SIA mechanism 
dominates. Developing a reliable experimental technique based on EDSR requires a tool 
that identifies the spin-orbit mechanisms contributing to spin-flip transitions and allows to 
establish, as applied to specific materials, the dominating mechanisms. To this end, we find 
the EDSR intensity for the Dresselhaus and Rashba models as a function of the magnetic 
field direction. Our results suggest that the angular dependence of the EDSR intensity is an 
unique characteristic of the various competing mechanisms of spin-orbit coupling contribut- 
ing to EDSR. 

In what follows we consider electrons confined in a parabolic QW along the z direction. 

Then the electron Hamiltonian is H = H + Hz + H S o, where 

h 2 k 2 mulz 2 , - 1 . 
= ~2^T + ~~t~ Z = 2^ <T9 ^ ^ 

are the orbital and Zeeman Hamiltonians, respectively. Here m is the electron effective mass, 

ujq is the characteristic frequency of the parabolic potential, k = — ?V + eA/hc, A is the 

vector-potential of the field B(8, if), and 9 and (p are the polar and the azimuthal angles of 

B. We have chosen parabolic confinement because it is known to be the only kind that can 



be solved exactly [12j for arbitrary B direction. The solution reveals the basic regularities 
of EDSR, including its dependence on the confinement strength. The spin-orbit interaction 
H so will be considered as a perturbation. 

Because Hq has quadratic form both in the momenta and coordinates it can be diago- 
nalized. Let us choose a new Cartesian frame where the z'-axis is parallel to B and the y' 
axis is in the QW plane, the Landau gauge A = (0,Bx',0), and introduce new variables: 
£ = x'cos7 — z' sin 7, and i] = x'sin7 + ,2' cos 7. Then H can be written as sum of two 



harmonic oscillators: 



H o = + «c)^C) ( 2 ) 

C=€,»? 



where = ^cos 2 7 + u 2 sm 2 (8 + 7) and = u^sm 2 ^ + u 2 cos 2 (8 + 7) are the 

frequencies, E{n^,n^) = S^^c( n C + 1/2) are the energy levels, and n^,, > 



121. Here 



We = eB/mc is the cyclotron frequency for B || z, and 7 is determined by the decoupling 
condition sin 27 = (uj q /uj c ) 2 sin [2(0 + 7)]. The operators and are defined by the 
following relations: C — Co = \Jh/2muj^(a^ + 6^) and A;^ = i^Jmuj^/2h(a^ — a^). The shifts in 
coordinates, £0 and rj , are related to the Landau momentum, k = k V ', as £0 cos 7 + ^0 sin 7 = 
X 2 k, where A = (hc/eB) 1 ^ 2 is the magnetic length. It is important that the operators 
of the kinetic momenta in the original frame, (k x ,k y ,k z ), which are used in the following 
calculations, can be expressed as linear combinations of the operators and a£. The 
coefficients depend only on the angles, 9 and ip, and the frequencies, and u v , and are 
independent of £ > V01 an d k. The frequencies of the coupled cyclotron-confinement modes, 
uj^(9), depend on the B direction. When oj c < uj 0i oj^{6) decreases from u c at 6 = to zero 
at 9 = 7r/2, while uj v (6) increases from too to {uj^ + uj 2 ) 1 ^ 2 . 

In presence of a high-frequency electric field, the term H so adds spin-orbit contributions 
to the time-dependent and time-independent parts of the total Hamiltonian H. The latter 
contribution leads to mixing of the spin levels and it is convenient to eliminate it, in the 
first order in H so , by a canonical transformation exp(T)P). The operator T is non-diagonal 
in the orbital quantum numbers (n^,n v ), and its matrix elements are 

n 5 , n a T n (l n„, a) = , — -, 3 

where a is the spin index. After the canonical transformation, the time-independent part of 
H conserves the electron spin projection on the magnetic field direction. 

Because the motion in the direction of the time- dependent electric field, E(t) || z, is 
confined by the parabolic potential, the time-dependent interaction ezE(t) is bounded. It 
does not depend on the spin, and the z-coordinate can be expressed in terms of and a v 



— \jfi/2muj£ sin(0 + 7)(ag + 



+^Jft/2muj ri cos(0 + "/)(a ri + a+). (4) 

The T transformation produces a commutator z so = [T,z]; hence, z acquires a spin- 
dependent part z so . Spin-flip transitions are induced only by the spin-orbit contribution 
ez so E(t) to the time-dependent part of H. 
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To find the intensities of the spin-flip transitions excited by an electric field applied in 
the z direction, one should calculate matrix elements of z so that are non-diagonal in spin 
projections and diagonal in the orbital quantum numbers. We will consider first the spin 
flip matrix element of z so in the case when spin-orbit interaction H so is dominated by the 
Rashba term 

Hr = a R (a x k y - a y k x ). (5) 

Of course, the parabolic confinement that is symmetric in z does not produce the Rashba 
term by itself. Therefore, we introduce it phenomenologically, e.g., as originating from the 
hexagonal symmetry of a wurtzite type crystal. A reliable estimate of the spin-orbit coupling 
constant can be obtained today only for the BIA mechanism that will be considered below. 
To simplify equations and elucidate the basic physics, we consider from now on the quantum 
limit, when only the lowest electron level = n v = is populated. 

A cumbersome algebra using several identities relating the frequencies u;| „ and the angles 
9 and 7 like cj| tan 7 = u>^ tan(0 + 7) results in a simple final equation for the matrix element 
of the spin-flip transition 

/t r I IN aR ^ c uj s (lj c - u s ) sin29 

The denominator of Eq. (0) can be rewritten explicitly as [ufao cos2 ^ ~ ^si^o + — u s))- 
The factor sin 29 that vanishes both for 9 = and 9 = it/ 2 reflects importance of a tilted 
magnetic field. 

The angular dependence of the EDSR intensity 1(9, (j)) oc |(f \z so \ 1)r\ 2 caused by the 
Rashba term is shown in Fig. la. Because of the poles of the denominator, the EDSR 
intensity increases when uo s approaches one of the eigenfrequencies; practically, for u c < uj 
only the pole uj^(9) = uo s is important. In the strong confinement regime, when u c , uj s ujq, 
uj£ pa uj c cos9 and becomes the cyclotron frequency of 2D electrons in a tilted magnetic 
field. The sharpness of the resonance peak is cut-off by a level width and also by the level 
anticrossing caused by the spin-orbit interaction. 

When the spin-orbit Hamiltonian H so is dominated by the bulk Dresselhaus spin-orbit 
interaction, the calculation of z so allows evaluating magnitudes of the EDSR for specific 
A3B5 compounds. In the principal crystal axes, the 3D Dresselhaus spin-orbit Hamiltonian 
Hd reads 

Hd = <^( cr ■ «0, where k x = k y k x k y — k z k x k z , (7) 



ky and k z can be derived by cyclic permutations, and 5 is a parameter. We have found a 
general expression for the matrix element (f \z so \ 1)d for a [0,0,1] QW in a A 3 B 5 crystal. 
The angular dependence of the EDSR intensity caused by the Dresselhaus term, 1(9, 0) oc 
|(T |^so| I)d| 2 ; is shown in Fig. lb. The expression for (f \z so \ 1)d simplifies in the strong 
confinement limit: 

£m cu c u; s sin g 

uj c — u s ) cos 6* sin2<^ — z(o; c cos 2 — u s ) cos 2y?] . (8) 

This equation describes a quasi-2D regime when IKd reduces to a 2D Dresselhaus term 
Hd = aD(a x k x — (Jyky) with a fl = —S(k 2 ) = —5mu /2h. 

Eq. (JHJ) has a pole at cj c cos 6* w ui s , similarly to Eq. (jSJ in the strong confinement regime. 
A distinctive feature of the BIA mechanism is a strong azimuthal dependence of the EDSR 
intensity, 1(9, (p), that possesses a four- fold axis symmetry. Remarkably, the contribution 
of the Dresselhaus term does not vanishes for an in-plane magnetic field, 9 = tt/2. In this 
geometry, the EDSR intensity does not depend on u s , which drops from Eq. (jHJ) because 
00^(71/2) = and shows an especially strong azimuthal dependence on the magnetic field 
direction, I(ir/2,(p) oc cos 2 2</>. 

Figure 1 shows a drastic difference in the EDSR angular dependences caused by the BIA 
and SIA mechanisms. The four-fold symmetry of the EDSR angular dependence should 
be broken when the contributions from both mechanisms to the EDSR amplitude are of 
comparable magnitude (This is also true for breaking the symmetry of the energy spectrum 
P. Fig. 1 shows convincingly that the angular dependence of EDSR is a powerful tool 
for identifying contributions of the different competing mechanisms of spin-orbit coupling. 

The efficiency of the BIA mechanism of EDSR can be evaluated by using the characteristic 
length, Id, that is equal to the matrix element of z so . Equation (jHJ) gives Id ~ 5m/2U 2 when 
all frequencies are of the same order of magnitude, u> ~ u c ~ uj s . We estimate Id ~ 10~ 9 
to 10~ 8 cm using a typical value m ~ 0.05mo for the mass and also S ~ 20eVA 3 for GaAs 
and 200 eVA 3 for InSb or GaSb Q|. It is much larger than the electron Compton length, 
— h/rriQC w 4 x 10 _11 cm, that plays the role of a characteristic length for EPR [14]. 
Therefore, electrical manipulation of electron spins is preferable to magnetic not only because 
it allows access to the spins at a nanometer scale but also because a larger coupling constant 
can be achieved. 
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However, there are several factors related to the electron confinement in a QW that 
can reduce Id- It is seen from Eq. (jHj) that the confinement frequency ujq appears in the 
denominator, therefore, Id includes a small factor uj c /uj <C 1. The factor uj c /uj reflects 
the deviation of the system from a strictly 2D geometry that is a critical condition for 
the gate- voltage controlled EDSR (the strict 2D limit corresponds to u — > oo). This 
factor was not really small in the Kato et al. experiment, uj c /ujq ~ 0.5, because a wide 
parabolic well with effective width about 50 nm and a strong magnetic field B = 6T were 
used jj]. In such a well ao ~ 0.3 x 10~ 10 eVcm which is much less than the typical value 
of the Rashba constant oir ~ 10" 9 eVcm for InAs based QWs and even larger values 
a^j ; ~ (3—6) x 10~ 9 eVcm were reported in Ref.[16]. This fact suggests that using asymmetric 
QWs should provide considerable advantages, and the corresponding length Ir may be larger 
;han Id- However, specific calculations of ccr depend strongly on the boundary conditions 



17j . and the dependence of oir on the QW width has not been investigated. 

Both Id and Ir are also reduced because of the spin-flip frequency uj s in the numerators 
of Eqs. (JBJ) and (JBJ). It introduces a numerical factor uj s /uj c = gm/2rriQ that is about 0.16 in 
GaSb and InAs and about 0.32 in InSb. Therefore, usually uj s /ujq rather than uj c /ujq is the 
factor controlling the intensity of EDSR. It originates because of the parabolic confinement 
in the z direction; a similar factor appears in the theory of EDSR for impurity centers jjj . 
For an in-plane magnetic field, motion in the B- direction becomes unrestricted, and that 
is why uj s cancels in Eq. (JHJ). In this case Id = (Sm/2h 2 )(uj c /uJo). Therefore, the orbital 
mechanisms of spin-orbit coupling can provide a strong EDSR only if the ratio u s /uj c is not 
too small. 

Another mechanism of EDSR, explored experimentally by McCombe et al. ^| for bulk 
InSb, is related to the anomalous coordinate f so = l^dc x k) introduced by Yafet Q|. The 
explicit expression for in the framework of the Kane model is Z so ~ Ti(\g\ / AnisEo) 1 ! 2 , where 
Eq is the forbidden gap It is rather large, Z so ^ 10~ 8 cm. However, the operator 

f so itself allows only the transitions at the combinational frequencies uj^ ± u s , hence, it can 
produce spin-flip transitions only in the second order of perturbation theory, in combination 
with some different mechanism, e.g., QW asymmetry. As a result, the corresponding length 
l r is much smaller than Z so . For a strongly asymmetric QW, the upper bound for this length 



is l T ~ (hu s / EG)L con f, where L con { = yh/mujQ is the confinement length. Unfortunately, we 
are not aware of any more specific experimental or numerical data on l r . 
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The spatial dependence of the electron g-tensor across the QW that allows efficient electric 
control of the Zeeman splitting jj, 0] also results in a spin-orbit coupling because Hz = 
fj,B{crg(z)B)/2 includes both the coordinate z and Pauli matrices. The magnitude of the 
electron spin coupling to the time dependent voltage V(t) is about ~ fiBV(dg/dV)B, where 
B = 6T. The derivative dg/dV is about dg/dV ~ g/Vo with g ~ 0.1 and Vq ~ IV, see 
Fig. 2B of Ref. I4. With V ~ Ew and w ~ 100 nm, the characteristic length l g caused by the 
spatial dependence of g is l g ~ XcgwB/Vo ~ 7 x 10~ 10 cm. It should be compared with Id 
that turned out to be negligibly small, Id ~ 10~ n cm, due to anomalously small ^-factor in 
GaAs (the ratio uo s /uo c = gm/2niQ < 10~ 2 ). Therefore, our estimates show that the gate- 
voltage manipulation of the electron spins achieved by Kato et al. was performed with a 
characteristic spin-orbit length I at the level of I ~ 7 x 10~ 10 cm. 

In most narrow gap semiconductor compounds with their typically large ^-factors, and 
especially in those with strong spin-orbit coupling, the BIA and SIA orbital mechanisms 
dominate the coupling of electron spins to a perpendicular electric field in moderate and 
strong magnetic fields, and the stength of this coupling is sufficient for electron spin ma- 
nipulation. Generally, however, relative contributions of the various spin-flip transition 
mechanisms depend strongly on the specific semiconductor material j2l| and the geometry 
of the QW. 

In conclusion, we have shown that the dynamic spin response to an electric field per- 
pendicular to the QW plane is controlled by the deviation of the confined electrons from 
strictly 2D behavior. Therefore, the response of the spin system to the gate voltage depends 
strongly on the ratio of the confinement layer thickness to the magnetic length, that should 
not be too small. Semiconductor compounds with large g-factors are highly advantageous for 
gate- voltage driven EDSR. Mixing the in-plane and perpendicular orbital motion is critical 
for EDSR, and for a (0,0,1) QW it requires that the magnetic field is tilted. The dependence 
of the spin-resonance intensity on the magnetic field direction with respect to the crystal 
axes is indicative of the role of the various mechanisms of spin-orbit coupling involved. The 
Dresselhaus and Rashba spin-orbit coupling mechanisms are efficient enough for electron 
spin manipulation by a time-dependent gate voltage. 
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FIG. 1: Angular dependence of the EDSR intensity 1(8, ip) for a (0,0,1) QW (in arbitrary units) 
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—0.17 (as in InAs) and oj c /ojq = 0.5 (as in Ref. |J|) were used for both mechanisms. To cut-off the 
pole in 1(6, ip) at uJi(0) = lo 2 s , an imaginary part iT with T = 0.08a; c was added to uj c . 
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